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Mode-mode coupling theory is presented for the resonant pumping via dynamical tunneling pro- 
cesses in a deformed microcavity. From the steady-state solution of the coupled differential equations 
of uncoupled chaotic modes and an uncoupled high-Q regular mode, pumping efficiency is obtained 
as a function of pump detuning, coupling constants and decay rates of the involved uncoupled modes. 
We show that the pump-excited chaotic modes as a whole can be regarded as a single pump mode 
with an effective decay rate and an effective coupling constant with the regular mode. We also show 
that the decay rate of the regular mode is enhanced by dynamical tunneling into all chaotic modes. 
Analysis method to obtain the effective coupling constant from the pumping efficiencies is presented 
for a two-dimensional deformed microcavity. 
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I. INTRODUCTION 

It is well known that a high-Q mode in a deformed 
microcavity can decay by chaos-assisted dynamical tun- 
neling. High-Q modes are usually localized in regular re- 
gions in a phase space surrounded by chaotic sea. With- 
out chaos-assisted dynamical tunneling the only way that 
the light in a high-Q mode can escape the cavity is by 
directly tunneling to the outside continuum. The output 
resulting from this is usually very weak. A strong di- 
rectional output can be obtained if there exists a chaos- 
assisted dynamical tunneling process in action: the light 
in a high-Q mode can undergo dynamical tunneling into 
nearby chaotic sea first and then go through chaotic ray 
dynamics until it escapes the cavity by ray refraction at 
particular positions and at particular angles. 

Quite recently, Yang et al. [l| has reported a resonant 
pumping experiment based on dynamical tunneling from 
chaotic sea to high-Q modes in a deformed microcavity 
laser. They introduced a collimated pump beam into the 
cavity by refraction in a time reversed way with respect 
to the directional output caused by dynamical tunnel- 
ing followed by chaotic ray dynamics. They observed 
that the pumping efficiency of a high-Q lasing mode at a 
longer wavelength than that of the pump is enhanced by 
two orders of magnitude whenever the pump is resonant 
with a high-Q cavity mode, which is localized in a regu- 
lar region in the phase space separated from the chaotic 
sea. Since the pump beam, injected by refraction, moves 
in the chaotic sea, the resonant enhancement must have 
come from the dynamical tunneling from the chaotic sea 
to the regular mode. 

In this paper, we present a mode-mode coupling the- 
ory for this pumping process. We first set up coupled 
differential equations between uncoupled chaotic modes 
and an uncoupled high-Q mode with the chaotic modes 
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driven by a pump beam in Sec.|TTJ From the steady-state 
solution of the coupled differential equations, we derive 
a formula for relative pumping efficiency as a function 
of pump detuning, coupling constants and decay rates of 
the involved modes in Sec. lIIIl We show that the pump- 
excited chaotic modes can be regarded as a single pump 
mode with an effective decay rate and an effective cou- 
pling constant g with the regular mode in Sec. lIII Cl We 
also show that the decay rate of the regular mode is en- 
hanced by dynamical tunneling into all chaotic modes. 
The same result is obtained by considering the eigenvalue 
problem associated with the coupled differential equa- 
tions in Sec.|lVl We finally describe analysis method to 
obtain the effective coupling constant g from the exper- 
imentally measured pumping efficiencies and apply it to 
the recent work by Yang et al. [l[ as an example in Scc.fVl 



II. MODE-MODE COUPLING MODEL 

Dynamical tunneling between an uncoupled regular 
mode and uncoupled chaotic modes can be modeled as 
mode-mode coupling, as depicted in Fig.[T] Here, the un- 
coupled modes (or states if you will) are not true eigen- 
modes of the system since the coupling between them 
are treated separately in our model - the term "mode" 
from now on actually means an uncoupled mode/state if 
not noted otherwise. Many chaotic modes (low-Q cavity 
modes) of decay rates 7„ (n = 1, 2, 3, • • • , N) are driven 
by a pump laser and these modes can then be coupled to 
a regular mode (a high-Q cavity mode) of a decay rate 
7 r by dynamical tunneling processes. 

Let us write the electric field of the nth chaotic mode 
as E n (x,t) — £„(£)/„ (x)e _lwt and that of the regular 
mode as E r (x,t) = £ r (t) f r (x)e~ %ut . Here £ n and £ r 
are their slowly-varying envelopes in time and /„ and f r 
are their normalized spatial mode functions, respectively. 
We assume that the chaotic modes have been obtained 
by a proper unitary transformation among them in such 
a way that they are orthogonal to each other. Hence, 
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FIG. 1: Uncoupled chaotic modes, excited by refractively in- 
jecting a beam of a pump laser, can couple to an uncoupled 
regular mode by dynamical tunneling processes. The nth un- 
coupled low-Q chaotic mode of an envelope £„ can be coupled 
to an uncoupled high-Q regular mode of an envelope £ r with 
a coupling constant g n (n = 1, 2, 3, . . . , N). The amplitude of 
a driving field is Eo with a coupling coefficient a n to the nth 
chaotic mode. 



there exists no coupling between any two of them. 



f*(x)f q (x)d 3 x = S pq where p, q = 1, 2, . . . , N 



The equations of motion for envelopes £ n and £ r are ba- 
sically those of driven coupled oscillators @: 



£n+ln£n — a n E —g n £ r , (1) 
N 

£ r + (j r + iA)£ r = y]g n £n, (2) 



where g n is the coupling constant (assumed real) between 
the nth chaotic mode and the regular mode. We as- 
sume that the chaotic modes are extremely lossy while 
the coupling of each to the regular mode is so weak that 
the oscillation between any chaotic mode and the regu- 
lar mode is over-damped: \g n \ <C j n - Effects of partial 
barriers [H are included in g n and j n . Symbol a n is the 
coupling coefficient of the external pump laser into the 
nth chaotic mode, depending on the position and the in- 
cident angle of the pump beam on the cavity boundary, 
and A = u — oj r , the detuning between the pump laser of 
frequency w and the regular mode of frequency u r . Since 
each chaotic mode is so broad in linewidth and since we 
are interested in a small frequency range (~ 7r) around 
the regular mode, we have neglected a detuning term for 
the chaotic mode in writing Eq. ([!]) . 



III. STEADY STATE SOLUTION 

By letting all g„ = 0, we simulate the non-resonant 
pumping. In this case in the steady state, we get 

£ n = C ^E ee £l 25 p (x) = J2 On« = £ p °./p(x), 

(3) 



where 



/pM = 4£0»(x) with t%=.[E\E°\* (4) 

C P n V n 



such that 



/;( X )/ P (x)d 3 x 



i. 



The quantity £® is the steady-state amplitude of the nth 
chaotic mode excited by the pump laser of Eq. Only a 
small portion of the chaotic modes that have appreciable 
a n/ln among a large number of chaotic modes make up 
the function / p (x). Therefore, / p (x) can be regarded 
as the (normalized) mode function of the nonresonant 
intracavity pump field or a pump mode in short. 

In the experiment of Yang et al. [lj , the position and 
the incident angle of the focused beam of the pump laser 
were adjusted in order to maximize the pumping effi- 
ciency This procedure is nothing but optimizing the set 
{«„} or the pump mode for the strongest tunneling to 
the regular mode. 

For nonzero g n , the steady-state solution is 
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where 
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7r + iA 
9ng 



7n(7r +iA) 



was neglected in Eq. © since 



iA^ 



o3« 



- £ 



9n 



(6) 

(7) 
(8) 
(9) 

(10) 



the magnitude of which is much smaller than and 
since \G\ is at most of order of unity as to be seen later, 
the quantity in the square brackets can be safely ne- 
glected. The pump mode is then obtained from 



£° 



(11) 
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where / p (x) is defined by Eq. (|4]). 

We are interested in the intensities, I r = \£ r \ 2 and 
I p = |£ p | 2 , for the evaluation of which we introduce the 
following quantities: 



and 



J_ = V — (— 

1 P ~ „ln\9 



(12) 



(13) 



We then obtain the intensities of the regular mode and 
the pump mode, respectively, as 



I r {5) = I 



G 



p 7r (1 + G) 2 + S 2 ■ 

7 (S ) - J° 1 + 62 
p[ > ~ p(1 + G) 2 +<5 2 ' 

G(2 + G) 



where 



I 



6 2 + (1 + Gf 



,~ P T. * = A/ 7r 



\£ 



(14) 



(15) 



(16) 



The relative pumping efficiency e with respect to the 
non-resonant pumping efficiency can be written as 



e(S) 



J [I p ffl|/ p (x)| 2 + I r (5)\fr(*)\ 2 ] I.A(x)| 2 d 3 x 
J mf p (x.)\*\Mx.)\*d*x 



(17) 

where /;(x) is the mode function of the lasing mode at 
which the pumping efficiency is measured in the work 
of Yang et al. [l| . We introduce mode-overlap factors f3 p 
and (3 r defined as 



|/ p (x)/ ; (x)| 2 ^x, p r = / |/ r (x)/ ; (x)| 2 ^x 



Using Eqs. (JT^J) and (|T5|) . we then obtain 

IrPp 



(18) 



(19) 



where a modified decay rate j' r , coupling efficiencies a, a' 
and a unity-peak-normalizcd lineshapc C(5) are defined 
as 



i r EE 7 r(l + G), 

G(2 + G) 
° - (1 + G) 2 ' 
G 



(1 + G)' 

(1 + G) 2 
<5 2 + (l + G) 2 



(20) 
(21) 

(22) 
(23) 



The coupling efficiency a is defined as the fraction of the 
intracavity pump power of non-resonant pumping that is 
reduced in the case of the resonant pumping, similarly 
defined to the case of tapered-fiber couplers Q. The 
other coupling efficiency a' is a fraction of the intracav- 
ity pump power of non-resonant pumping that is trans- 
ferred to the regular mode in the resonant pumping. The 
physical meaning of "f' r is to be discussed in Sec. MI Cl 



A. Effect of linewidth of the pump laser 

The above results were obtained for a monochromatic 
pump laser. If the pump has a linewidth jl comparable 
to or larger than the linewidth 7 r of the high-Q cavity, 
we should average the results over the spectral lineshape 
of the pump 



C l {uj') 



7i 



(a/ - w) 2 + 7 2 



with oj interpreted as its center frequency. The lineshapc 
function in Eq. (|23[) then becomes 



C(5) tt(1+G) 



(1 + G)/tt 



(7i/7r)/7T 



dS', 



5' 2 + (1 + G) 2 (6' - 5f + ( 7i /7r) 2 

(24) 

where 6' = (u/ — uj r )/"/ r and 5 = (oj — u r )/^ r . 

If 7l <C 7 r , the second integrand can be approximated 
as a delta function and thus we recover Eq. (|23|) . If 7l 3> 
7r = 7r(l + C) on the other hand, the first integrand is 
like a delta function. 



C(S) — » tt(1 + G) 



1 



(7i/7r) 



7T <5 2 + (7i/7r) 2 



(25) 



The linewidth of the resonance appears to be that of the 
broad pump laser and there is also a reduction factor 
l' r /lL *C 1. Eq. (|19[) is now replaced in this limit with 



e(S) 



1^ 

1L 



aC L (6) 



+ 2p^ a 'C L (5). (26) 



B. Multiple interference picture 

The steady-state solution given by Eqs. © and © can 
also be obtained from consideration of multiple interfer- 
ence of coupled fields. Without the coupling, g n = 0, the 
chaotic modes are excited to £ n = £^ with no excitation 
of the regular mode, £ r = 0. In the presence of the cou- 
pling, g n ^ 0, the chaotic modes can tunnel to the regular 
mode. Let us consider the zero detuning case, A = 0, for 
simplicity. Equation ^ dictates that the amplitude of 
the regular mode tunneled from the nth chaotic mode is 
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with g n /j n as a tunneling coefficient composed of the 
coupling constant g n and the decay rate of the destina- 
tion mode 7 r . Summing up the contributions from all 
chaotic modes, we obtain the first round contribution to 
the regular mode: 



4 1] = E ?e 

Jr 



(27) 



This amplitude can tunnel back to the nith chaotic 
mode, the amplitude of which is given by 

\ 1™, J 

where the tunneling coefficient (—g m /"fm) is dictated by 
Eq. fl]). Note the decay rate in the denominator is 7 m 
since the destination this time is the mth chaotic mode. 
A fraction of this amplitude given by 



9m \ c(l) _ 9m c(l) 
I c m c r i 

7r / 7m7r 



can tunnel back to the regular mode. Summing up the 
contributions from all chaotic modes, we obtain the sec- 
ond round contribution: 



£^ 



7r7ri 



Extending this argument, we easily find the fcth round 
contribution to be 



°r 

£i k) 



9n p(k) 
In 



(28) 



Summing up all multiple round contributions, we then 
obtain the amplitude of the regular mode in the form of 
multiple interference. 



£ f — E^. ^ ~\~ ^ ~\~ £!^ ^ ~~\~~ 
= £^(l-G + G 2 ) 



& _E„Sn^/7r 



1 + G 



1 + G 



recovering Eq. with A = 0. Likewise, the amplitude of 
nth chaotic mode is obtained by summing up all multiple 
round contributions in the form of multiple interference: 

£ - £° + £ {1) + £ {2) + ■ ■ ■ 



- £ u - 



In 

1, 



(£^+£^+£^ + ---) 



cO _ 9n c 



which is nothing but Eq. ^ with A = 0. 



C. Physical meaning of g, 7 P , G and 7,'. 



In Eqs. ([7]) and (fT2|) , g and -f p were defined, respec- 
tively, without identification of their physical meanings. 
We show below that g is the coupling constant between 
the regular mode of £ r and the pump mode of £ p and 
that 7 P is the decay rate of the pump mode. 

Consider the first round contribution to the regular 
mode in Eq. (f2"7) . By using the definition of g of Eq. J7]), 
this can be rewritten as 



1, 



■£° 



(29) 



indicating that it can be regarded as a result of tunneling 
from the pump mode of £ p with a coupling constant g. 

We can compose the first round correction in the pump 
mode from the first round correction in the nth chaotic 
mode. From Eq. © with A = 



£° 
1 + G 



£n ~ G£® + G 2 £® 



and thus 



* -GE2 



The first round correction in the pump mode is then given 



as 



£«(x) = ^£W/ B (x)~-G^^(x) 

n n 

= -G£°f p (x) 
From Eqs. (JT2J) and {jT3j) we have 

,2 



y_ = y al = lrG 

lp „ In 



(30) 



(31) 



Therefore, 



£P = -G£° = 



Iplr 



9 " 



9_ F (i) 

°r 



indicating that the first round correction £^ in the pump 
mode can be regarded as a result of tunneling from the 
regular mode of £^ with the same coupling constant g 
and a decay rate j p associated with the destination mode 
or the pump mode. 

Our identification of g and 7 P is not restricted to the 
first round corrections. By using the physical meaning 
of g and 7 P , the second round correction to the regular 
mode would be given by 

c(2) _ 9_c(l) 

c r - c p , 

Jr 



which can be simplified as 



£M 



9_ 



9_ 

1p 



5 



which is consistent with Eq. (|29[) obtained without relying 
on the present identification of g and j p . Likewise, we 
expect the second round correction to the pump mode be 
given by 



p(2) _ 
°p - 

which can be simplified as 



Irlp 



Ip 



This consideration can be extended to all higher-round 
correction terms. Summing up all correction terms, we 
obtain 



f 



c p 
1 + G ' 



which is the same as Eq. , also obtained without using 
the present identification of g and 7 p . Therefore, we con- 
clude that g is the coupling constant between the pump 
mode and the regular mode and j p is the decay rate of 
the pump mode. 

The physical meaning of G and 7^ can be obtained as 
follows. According to the cavity quantum electrodynam- 
ics (QED), when a high-Q oscillator with a decay rate 
jh is coupled with a coupling constant g to a low-Q os- 
cillator with a decay rate jl (3> Jh), total decay rate of 
the high-Q mode is enhanced as @ 



lH = 1H + 



£_ 

1L 



1H 1 



1H1L 



) 



(32) 



This phenomenon is known as the enhanced spontaneous 
emission in the cavity QED. Comparing Eqs. (|3"2"j) and 
(|3"Tj) , wc can identify that the righthand side of Eq. (f3"Tj) 
is nothing but the increment in the decay rate of the 
regular mode induced by tunneling into all chaotic modes 
and that G is just an enhancement factor with respect 
to the uncoupled decay rate j r . Total decay rate of the 
regular mode modified by the tunneling is then given by 
7 r (l + G), which was defined as 7^ in Eq. (|20|). In fact, 
this total decay rate is correctly reflected in the lineshapc 
C{5) of the resonant pumping efficiency in Eq. ()23l) . 

Equation (|3Tj) shows that the total tunneling rate of 
the regular mode into all chaotic modes is equal to g 2 /j p , 
which is an effective tunneling rate of the regular mode 
into the pump mode. This identification is in fact con- 
sistent with the Fermi golden rule : 

W=^-\{p\V\r)\ 2 6 E , 

where r and p stand for the regular and the pump modes, 
respectively, V is the interaction Hamiltonian, and 5e is 
the density of states. If V is known, we can calculate its 
matrix element. Let us denote it as 




pump 



FIG. 2: The regular mode can couple to all chaotic modes 
(n = 1,2, ■•■ , N) with an individual tunneling rate gn/jn- 
We can view this by introducing a pump mode, a specific 
collection of chaotic modes excited by the pump laser, which 
has an effective coupling g with the regular mode and an 
effective decay rate y p . 



The density of states is obtained by counting the number 
of modes per unit energy interval or by taking the inverse 
of the energy interval associated with the pump mode 



7 P 2 



((J - UJ r ) 2 + il 



1 



27rfi,7 P 



The extra 1/2 factor comes from the fact that only one 
polarization direction is allowed for the tunneling process 
out of two possible polarization directions. Therefore, the 
tunneling rate into the pump mode given by the Fermi 
golden rule is 

\g\ 2 



h ^ ^ 2nhrj p 7 p 



Lastly, with the above identification of 7 P and 7^, the 
meaning of the factor 7 P /7^ appearing in the pumping 
efficiency formula, Eq. (|19p. becomes clear. It is an 
intensity build-up factor. In non-resonant pumping, the 
intensity of the pump mode is inversely proportional to 
its loss rate j p , whereas in the resonant pumping the 
intensity of the regular mode, as a true eigenmode of the 
system, is inversely proportional to its total decay rate 
Y r . So, the factor j p /j' r measures how much intensity 
build-up is enhanced by the resonant pumping compared 
to the nonresonant pumping. 



(j>\V\r) = kg 
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D. Validity of rate equation model 

One may wonder whether the pumping efficiency for- 
mula of Eq. (fT9|) might also be obtained from rate equa- 
tions for intensities not electric fields. In order to address 
this question, let us consider a two-mode problem for sim- 
plicity, in which the intensities I r and 7 C in the regular 
and chaotic modes, respectively, satisfy the following rate 
equations. 



Ic 

i, 



R - 2 7c 7 c 



-2 7r 7 r 



7r 



2ff 2 



7r 



7c 



(33) 
(34) 



where R is a pumping rate for the chaotic mode. The first 
term in the square bracket is the energy transfer from the 
chaotic mode to the regular mode while the second term 
is that from the regular mode to the chaotic mode. Since 
G = g 2 1 {idr), we can rewrite these equations as 

i e = R - 2 7c (l + G)I C + 2 lr GI r , 
i r = -2 7r (l + G)I r + 2 7c G7 c . 

Since rate equations are valid in the limit of large decays, 
the above equations are valid only when G<1. 

Nonetheless, the steady-state (7 C = I r = 0) solution is 
obtained as 



Ic 
Ir 



= I 



l + G 



I 



1 

o7c 



2G 
G 



7r 1 + 2G 



7°(1-G), 



7 C °^G, 

7r 



(35) 
(36) 



where 7° = 77/(2 7c ). Unlike Eqs. (TJU) and ([15]). the re- 
sults of the mode-mode coupling model, the rate equation 
model shows that the power removed from the chaotic 
mode would be all transferred to the regular mode. More- 
over, from Eqs. ([33)) and (|34|) we obtain a statement of 
energy balance: 



R 

2 7c 7 c + 2 7r 7 r = 2 7c 7° 
constant = (energy loss) 



(pumping) 



which is consistent with Eqs. (|35|) and (|36|) . 

This energy balance is not satisfied in the mode-mode 
coupling problem of Eqs. (jTJ and @, which is simplified 
for two modes of £ c and £ r on resonance as 

£ c = a c E - 7c £ c - g£ r , 

£ r 7 r^r ~l~ g£c j 



which can be rewritten in terms of intensities as 



7 C = 2U[a c E £*} - 2 7c 7 c 
i r = -2-i r I r + 2g%t[£* r £ c 



2g^[£* r £ c 



(37) 
(38) 



where . .] represents a real component. Comparing 
these equations with Eqs. ([3"3"f and (|3~4"]) . we identify the 
following correspondence: 



2ft[a c E £*} o R, 

2a 2 2a 2 
2gW£* r £ c \ f> J-I c -J-I r 

Ir Ic 



(39) 
(40) 



However, the pumping term 2^t\a c E^£ *} is not constant 
but dependent on £ c . This dependence is caused by the 
interference of the pump field with the field of the chaotic 
mode. Using Eq. ^ , we can rewrite the pumping term 
as 



2R[a c E £*} = 



2M[a c E £°} _ R 



l + G 



l + G ' 



(41) 



which decreases as G increases. In fact, from Eqs. (TT4")) 
and (|15[) in the steady state, we have 



(pumping) 



2^[a c E £* c ] 
2 7c 7 c + 2 7r 7 r 



2 7c 7 r ° 
l + G 



from which we identify Rq = 2j c I®. 

The rate equation model and the mode-mode coupling 
model should give the same result in the limit of large 
decays or G C 1. However, they give seemingly differ- 
ent results, Eqs. ([33]) and (|3"6"j) compared to Eqs. ([T4"|) and 
([151) . The difference originates from the assumption that 
the pumping rate would be constant in the rate equa- 
tion model. We have found in Eq. (|4"Tj) that it varies as 
(1 + G) . In fact, if we replace the pumping rate R in 
the rate equation model with R = 2 7c 7" /(l + G), we then 
obtain 



Ic 
Ir 



I°c 
I 



1 



1 + 2G 

o7c 



:7 C °(1-2G), 



7r (1 + G)(l + 2G) 



7^G , 

7r 



(42) 
(43) 



which is now perfectly consistent with Eqs. (fT4")) and (|15l) . 
the results of the mode-mode coupling model, in the limit 
of G < 1. 



IV. EIGENVALUE SOLUTION 

The decay rate 7r is the decay rate of the regular mode 
when it is uncoupled to other modes. However, the de- 
cay rate directly measured in experiments is not 7r , but 
a modified decay rate due to the tunneling into all 
chaotic modes. By physical considerations in the preced- 
ing section, we identified 7 £ defined in Eq. (|20|) is this 
modified decay rate appearing in the lineshape of the 
regular mode. 

The expression for the modified decay rate or a total 
decay rate of the regular mode as a true eigenmode of 
the system can also be rigorously derived by solving the 
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mode 
I 


e 

(exp) 


/ 

7r 

(exp, in 10 9 s _1 ) 


(sim) 


7r 

(in 10 9 s" 1 ) 


G 


9 

(in 10 10 s- 1 ) 


a 


/ 

a 


7rG 

(in 10 9 s- 1 ) 


1 


25 ±4 


2.0 ±0.6 


1.00 ±0.05 


1.7 ±0.5 


0.19 ±0.09 


2.3 ±0.7 


0.29 ±0.11 


0.16 ±0.06 


0.32 ±0.18 


2 


60 ±8 


7.9 ±2.3 


0.9 ±0.1 


4.5 ± 1.9 


0.75 ±0.53 


7.6 ±2.9 


0.67 ±0.20 


0.43 ±0.17 


3.4 ±2.6 


3 


22 ±4 


20 ±6 


0.6 ±0.1 


15 ±5 


0.30 ±0.17 


8.7 ±2.9 


0.41 ±0.16 


0.23 ±0.10 


4.5 ±2.9 


4 


5.5 ±0.7 


260 ± 80 


0.4 ±0.1 


190 ± 70 


0.36 ±0.27 


34 ± 14 


0.46 ±0.21 


0.26 ±0.14 


68 ±57 


5 


2.7 ±0.4 


790 ± 230 


0.3 ±0.1 


440 ± 190 


0.80 ±0.62 


77 ±34 


0.69 ±0.21 


0.44 ±0.19 


350 ± 310 



TABLE I: Decay rate y r of the uncoupled regular mode, the tunneling-induced enhancement factor G in the modified decay 
rate ~/' r of the regular mode, the effective coupling constant g for the pump mode, the coupling efficiencies a and a', and the 
total tunneling rate 7 r G are all obtained from our analysis. 



eigenvalue problem of Eqs. ([T]) and <j2j) on resonance (A 
0), which can be rewritten in a matrix form as 



dt 



E = -rE ± E Q A , 



(44) 



where 
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-.91 • ■ 


• —9N 







E = 


£i 


, r = 


9i 


7i •• 


• 


, A = 


Ol 




£n 




gN 


•• 


• 7-/V 




a N 



The steady-state solution in Eqs. ([5]) and (J6)) are obtained 
by letting the lefthand side vanish and multiplying the 
inverse matrix r -1 to both sides: 



E = £ r _1 A 



(45) 



For transient response of the system we also need to con- 
sider the homogeneous equation E = — TE. The solution 
of the homogeneous equation determines the spectrum of 
the system. For a trial solution, we assume E cx e~ At , and 
the result, TE = AE, is nothing but an eigenvalue equa- 
tion with the eigenvalues obtained from dct(r — AI) = 0. 
The determinant can be arranged as 



det(r - AI) 



N 



J(ln A) 



N 



(in - A) 



(46) 



So, det(r - AI) = if A = 7n (n = 1, 2, . . . , N) or if the 
quantity in the second bracket vanishes. Since 7„ 3> A ~ 
7 r , the expression in the second bracket is simplified as 



N 



9l 



A~ 7r +^^= 7r (l±G). 



(47) 



The eigenvalue A = 7„ has the eigenvector Ej = 
5j„ = ip„ (there are N such eigenvectors), which is the 
original uncoupled nth chaotic mode. The eigenvalue 
A = 7£ = 7, (1 ± G) corresponds to the modified decay 
rate of the regular mode. This result is consistent with 
the one obtained from the cavity-QED and the lineshapc 



considerations in the preceding section. The eigenvec- 
tor ip' r corresponding to this eigenvalue is the modified 
regular mode and it is easily obtained as 



ip' T = V'r 



(48) 



up to the first order of g n /ln- Here ij) r is the uncoupled 
regular mode. The modified regular mode, which is what 
we measure in the spectrum as a true eigenmode of the 
system, has a small contribution from each chaotic mode 
with a relative amplitude of (g n /jn) *C 1. 



V. ANALYSIS 

In the experiment by Yang et al. [l[ , the relative pump- 
ing efficiency was measured on resonance, the formulae 
for which are 

e(0) = [1 - a(G)] + ^a'(G), for lL « j' r 



% 
1L 



a(G) 



^«'(G),for 7 L»7; 



The pumping efficiency e(0) depends on the following pa- 
rameters: 7 r , 7r, 7 P , 7l, G, P p and /? r , among which j p , j3 p 
and f3 r are obtained from numerical simulations (see be- 
low) and 7^ the pump laser lincwidth and 7^ the observed 
lincwidth of the regular mode are obtained from experi- 
ment. There are left two unknown parameters j r and G. 
There is one more relation to use, that is, Y r = j r (l + G). 
Therefore, we can simultaneously solve two equations to 
obtain j r and G. Once they are found, we can then ob- 
tain the effective coupling efficiency g = j p J2 n 9nlln- 

The decay rates of the pump mode 7 P can be calculated 
from the following ray simulation Q : a bundle of rays is 
initially prepared in the pumping position with a pump- 
ing angle which correspond to the experimental pumping 
condition and the path length L p of the bundle of rays 
before it escapes the cavity is calculated. The decay rate 
7 p is then simply c/ (2mL p ) with c the speed of light and 
m the refractive index of the cavity medium. The overlap 
factor p p between the pump mode and the lasing mode 
of interest can also be obtained from a similar ray simu- 
lation. The overlap factor (3 r between the regular mode 
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and the lasing mode can be obtained from wave calcula- 
tion, by using polar-angle-averaged mode distributions. 
Table 1 summarizes the results of the present analysis 
applied to the experiment of Yang et al. [l| . 

VI. CONCLUSION 

We have developed a mode-mode coupling theory for 
the resonant pumping via dynamical tunneling in a de- 
formed microcavity laser. We derived a formula for the 
pumping efficiency, which is readily measured in exper- 
iments, as a function of pump detuning, coupling con- 
stants between an uncoupled high-Q regular mode and 
uncoupled chaotic modes corresponding to the chaotic 
sea surrounding the regular region supporting the reg- 
ular mode, and decay rates of the regular and chaotic 



modes. Analytic expressions for coupling efficiency of 
the pump into the regular mode is derived. It is shown 
that the pump-excited chaotic modes collectively couple 
with the regular mode with an effective coupling constant 
and a single effective decay rate. We show that the cou- 
pling efficiency and the effective coupling constant can be 
obtained from the observed pumping efficiency on reso- 
nance by using our theory. We applied our theory to the 
experiment by Yang et al. and obtained the coupling ef- 
ficiency as well as the effective coupling constant as an 
example. 
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